arXiv:1505.07320vl [physics.optics] 27 May 2015 


Optical vortices: the concept of topological potential and analogies with 

two-dimensional electrostatics 
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We show how the phase profile of a distribution of topological charges (TC) of an optical vortex 
(OV) can be described by a potential analogous to the Coulomb’s potential for a distribution of 
electric charges in two-dimensional electrostatics. From what we call the Topological Potential (TP), 
the properties of TC multipoles and a 2D radial distribution were analyzed. The TC multipoles have 
a transverse profile that is topologically stable under propagation and may be exploited in optical 
communications; on the other hand, the 2D distributions can be used to tune the transverse forces 
in optical tweezers. Considering the analogies with the electrostatics formalism, it is also expected 
that the TP allows the tailoring of OV for specific applications. 

PACS numbers: 42.50.Tx, 42.25.-p 


Vortices are ubiquitous iu uature, aud may appear 
wheu a physical field rotates arouud au axis. They play 
au importaut role iu fluid turbuleuce [1], atmospheric 
pheuomeua [2], supercouductors/superfluids mm aud 
optics [5]. ludeed, optical vortices (OV) are uuique be¬ 
cause the electric field of light evolves liuearly iu free- 
space, aud therefore, OV form a prototype system for uu- 
derstaudiug the geueral properties of vortices [6]. Siuce 
the semiual work by Alleu et al [7] OV are subject of 
iuteuse research, beiug applied iu areas as diverse as opti¬ 
cal tweezers, laser traps aud atom guides HHini , to excite 
surface plasmous IU, aud iu classical aud quautum com- 
muuicatious mils]. 

A quautity preseut iu OV is the topological charge 
(TC), which measures the uumber of twists the electro- 
maguetic field gives arouud a giveu axis iu a waveleugth. 
The TC is a uou-zero iuteger for a vortex, aud character¬ 
izes their uou-trivial topology m However, while the 
TC determiues the vortex homotopy class, it does uot 
determiue the geometry; OV cau possess the same TC 
but distiuct geometries m- 

Iu this work we demoustrate that a distributiou of TC 
iu OV cau be described iu a similar way to that of a two- 
dimeusioual (2D) distributiou of electric charges. The 
formalism hereiu preseuted may be used to desigu aud 
to uuderstaud paraxial OV iu uuusual geometries aud 
also should give iusights for vortex properties iu other 
physical systems. 

Cyliudrical coordiuates, where r = r(r, 0, 2 ;), are cou- 
sidered throughout the paper. Iu OV, the Poyutiug vec¬ 
tor (aud also the beam wavefrout) rotates arouud the 
propagatiou axis r = 0. The phase of the electric field 
E (r) is uot defiued at r = 0, so E (r = 0) = 0. Wheu 
au optical vortex is au eigeustate of the light orbital au- 
gular momeutum (0AM), the TC is equal to the 0AM 
per photou m A poiut TC of charge q is associated 
with a cyliudrical vortex profile, aud appears iu E (r) as 
a phase term where (p is the azimuthal augle. The 
aualogy betweeu au electrostatic charge distributiou aud 
a TC arraugemeut is cousidered here by uoticiug that 


addiug two TCs qi aud q 2 at the same poiut origiuates 
a phase e*Oi+92)0 gQ^ compositiou of TC is addi¬ 
tive iu the correspoudiug electric field. Also, siuce the 
total TC is usually couserved iu physical processes, OV 
are resilieut to atmospheric turbuleuce m aud decoher- 
euce m- Aualogously, the effect of two electric charges 
placed at a poiut is also additive, aud the total electric 
charge is a couserved quautity. 

The azimuthal phase term coustraius the beam iuteu- 
sity profile uear r = 0 at ^ = 0. We assume that the 
vortex is iuside a smooth euvelope A (r) arouud the ori- 


giu such that V^A(r) 


r=0 


0 , where Vx is the gradi- 


eut takeu over the plaue r<p. The electric field E (r) for a 
mouochromatic aud homogeueous liuearly polarized light 
beam propagatiug aloug z axis ou the paraxial regime is 
represeuted by 


E (r) = x£ (r) = xf (r) (r) 

( 1 ) 

and must satisfy the slowly varying envelope approxima¬ 
tion (SVEA) 


(r) = (r). (2) 

Equatiou 0 implies that for auy q dX z = {), where 
§^£ (r) = 0, Vi (/e*«‘^) = 0. So, / = rl«l at r « 0. 
For g > 0 , one may verify that the vortex term 
is the polar represeutatiou of the complex uumber 
while for q < 0, = (rt*) Therefore, we 

assume that a TC of q > 0 cau be represeuted by u^. 
However, siuce the topological properties are iuvariaut 
uuder small coutiuuous geometric perturbatious, aud 
{u — 6 uY will roughly represeut the same vortex for a suf- 
ficieutly small Su. Also, if q = qi q 2 , we may cousider 
that ^ {u — {u — 6 u 2 Y^ for small 6 ui {i = 1 , 2 ) 

ca. A possible geueralizatiou of these trausformatious 
at z = 0 is giveu by 


N 

exp V (r, 0, 0) = ]^ (44 — = i?exp {i^r) • (3) 

i=l 



2 


In the limit g ^ 0, ^ oo with q/N constant, it 

is obtained that V (r, (f>) = f da'p {u') log {u — u'), where 
p{u') is an effective TC density per unit area at u' = 
r'e^^ This expression for V is formally identical to the 
2D electrostatic potential [iniiini- By extending the pre¬ 
vious calculation to include g < 0 , we have the more 
general expression, 

y(r, 0, 0) = J da' [\p{u')\\og\u — u'\ ip{u') 8iTg{u — u')] . 

(4) 

Considering the similarity with the 2D electrostatics 
Coulomb potential, we interpret V (r, 0, 0) as a Topolog¬ 
ical Potential (TP) due to an arbitrary TC distribution 
over the plane z = 0. Accordingly to this interpretation, 
the TP of a TC distribution on a plane is related to the 
potential associated with infinite linear distributions of 
electric charges. Of course there are some important dif¬ 
ferences. For example, the fact that R = exp {fHe [V (r)]} 
must be finite everywhere implies that [V (r)] must 
be always positive and cannot include the TC signal. In 
other terms, the intensity profile is not sufficient to de¬ 
termine the topological properties of an optical vortex. 
However, the phase profile, 3m [V (r)], contains all the 
topological properties, and it also constrains the inten¬ 
sity profile. It must be remarked that although fHe [V (r)] 
is only an approximation valid near r = 0, Jm [V (r)] 
is exact. Another important difference between the TP 
and electrostatics is that the TC is necessarily a discrete 
quantity and, as will be shown later, the continuum gen¬ 
eralization implied by Eq. @ indicates that in the gen¬ 
eral case p {u') is not a direct map of the TC distribution. 

As a remark to our interpretation of Eq. (^, the reader 
should notice that it was insinuated in m that the vor¬ 
tex phase could be understood as a potential, and the TP 
introduced in this work extends this concept by adding 
the spatial structure. 

Another important reason to consider Eq. @ as a 
potential comes from the expression for the transverse 
paraxial momentum density of light, pj_. It is possible 
to write in our notation that [7] 


which by using the Cauchy-Riemann conditions gives 


' core 


( 9 ^') = ^ y da'\p{u')\avg{u-u') 


(6) 


e(0) 


This sort of relation between the phase profile and vor¬ 
tex core geometry was empirically inferred from experi¬ 
mental data [53] and verified under more general condi¬ 
tions in [24|. However, in the present formalism it arises 
naturally. Since we did not consider aperture effects, Eq. 

is expected to be valid when p is concentrated near 
r ~ 0. If p (r) has the same sign for all r, Eq. ^ simpli¬ 
fies to Vqqyq = |9^t/ 90|, where the azimuthal derivative 
is called Local Circulation (LC) and has an important 
meaning. It is locally proportional to the classical 0AM 
density and for a point TC it gives exactly the total 0AM 
per photon (TC) g [16]. The LC in the present context 
is the local 0AM as discussed in m- The vortex radius 
may be adjusted by tuning the LC to obtain OV with de¬ 
signed shapes [23l|24]. It is also possible, by distributing 
the TCs over given geometrical patterns to obtain more 
general OV profiles as lines, corners and triangles m- 
These procedures can shape the OV profile at 2 : = 0. 
However, the necessary conditions for beam profile sta¬ 
bility under propagation for beams produced via a TCs 
distribution remains as an open question. The connec¬ 
tions between the TP and the theory of spiral light beams 
[25] may provide an answer to this point. But this dis¬ 
cussion is outside the scope of the present paper. 

An electrostatics-related feature from Eq. Q is the 
existence of TCs multipoles. Since multipoles are usu¬ 
ally composed of oppositely charged TCs, they may an¬ 
nihilate under propagation and may be unstable [26II28] . 
However, by expanding the integration kernels of Eq. 0 
in terms of circular harmonics, one may verify that a vor¬ 
tex may be represented in general as V (r, 0, 0) = Aq (r) + 
iBo (r, (j)) + [Aj (r) cos + Bj (r) sin (j^)], where 

Aj (r), Bj (r) are complex functions. We remark that 
point multipoles with Aj and Bj proportional to 
are not meaningful for OV because these terms imply 
that when r ^ 0, both fHe [V (r)] and 3m [V (r)] diverge. 
We consider that a pure TC multipole is represented by 


Px(x3m{VxE(r)}/(r). (5) 


= asm{j(f)^ (3)/j, (7) 


In a semiclassical interpretation of Eq. <§ , we may 
understand the intensity profile / (r) as the probability 
of finding a photon at a position r, while 3m{V±V (r)} 
represents the local transverse momentum of the photon. 
Thus, considering that such light beam can transfer this 
momentum to an object, the associated force term would 
be proportional to the gradient of the TP. 

An application for Eq. 0 is to shape the core of an 
optical vortex at 2 ; = 0. Assuming that at r ^ 0 the 
vortex envelope is Gaussian pg, ^(r) - e / 2 , the dark 
core profile of the vortex is described by r • (r) = 0 , 


where a is constant and j 7 ^ 0. Eor an azimuthally peri¬ 
odic solution of Eq. 0 j must be an integer. A fractional 
j leads to a line of phase discontinuity similar to those 
of [29l [30] . /3 is an orientation offset while a determines 
vortex core local radius via Eq. (|^, such that 

rcore (</>) = Vlacos (jc/) +/3)|. ( 8 ) 

To produce and characterize some experimental con¬ 
sequences of the previous discussion, we modulated the 
wavefront of a 800 nm fiber-coupled laser diode by using 
a liquid crystal phase-only spatial light modulator (SLM) 
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in one arm of a Michelson interferometer, with the same 
experimental setup and detection scheme as previously 
described in m- Unless otherwise stated, the data was 
collected at the SLM image plane (^ = 0 cm). The phase 
profiles on the SLM were composed of a carrier wave, a 
circular apodization with a fixed radius and the phase of 
interest. 

We show in Figs, [^(a, d, g) the experimental intensity 
profiles for TC multipoles of order j =1, 2 and 3 by ap¬ 
plying the phase profile of Eq. 0 to the SLM. The solid 
lines represent the expected core profile as given by the 
LC. Blue and red lines surround, respectively, regions of 
negative and positive LC. The solid lines have only the 
maximum radius as an adjustable parameter, and since 
(a = 40 for all j, the same value was used for all curves. 
Figs. (b, e, h) show the LC, as determined 

from the experimental data m, and Figs, [^(c, f, i) ex¬ 
hibits the expected LC according to Eq. Q. A good 
agreement is found between the theoretically expected 
results and the experimental findings. The disagreement 
exists only at the darkest regions near the profile center, 
where we were not able to properly measure the phase. A 
technical aspect which is worth noticing is that the inten¬ 
sity profile of multipoles is very sensitive to the spatial 
filter iris transverse position, and misalignments makes 
the lobes profile nonsymmetrical. 



Figure 1. Multipoles of TC at z = 0 cm. The rows correspond 
to the data for a dipole (a, b, c), quadrupole (d, e, f) and an 
hexapole (g, h, i). In columns we display the beam amplitude 
profiles (a,d,g), the experimental LC (b, e, h) and theoretical 
LC obtained from Eq. 0 (c, f, i). The solid lines in (a, d, 
g) corresponds to the expected OV core profile from Eq. 
and their colors (red, blue) represent the enclosed TC sign 

(+5 ■)• 

An important property of TC multipoles is that their 
LC is stable under propagation, as can be seen by vary¬ 
ing the position of the CCD along the ^ axis. Exper¬ 
imental amplitude and LC at different values of 2 : are 


shown in Fig. for j = 4. It may be observed in the 
amplitude profiles. Figs. |^(a-d), that pairs of amplitude 
lobes with opposed LC signs annihilate under propaga¬ 
tion. The resulting bright spots are located at zero LC 
regions. Creation and annihilation of oppositely charged 
OV pairs under propagation are well established in lit¬ 
erature mi nil ED EH, but to our knowledge the pre¬ 
vious descriptions were always associated with, respec¬ 
tively, creation and destruction of TC. In the case we de¬ 
scribe here, the beam’s topological structure is preserved 
under propagation, as can be seen from the LC in Fig. 
(e-h). For negative z, the amplitude lobes rotate in the 
opposite direction to that shown in Figs. i(a-d). 


z=0 cm z=2 cm z=4 cm z=6 cm 



Figure 2. Propagation of a vortex octupole {j = 4) in free- 
space. In the top row (a-d) it is displayed the amplitude 
profile at increasing propagation distance from the SLM image 
plane (z = 0 cm), while in the bottom row (e-h) is shown the 
experimental LC. Notice that the LC remains stable under 
propagation. 


Another important property of the TP, Eq. 0, is 
that it satisfies a Gauss law inside a region D for the 
total enclosed TC, 


(£ dl • VU (r) = 27ri f da'p (r'). (9) 

JdD J D 

Equation (|^ is a generalization of the usual expression 
for the winding number over the phase profile, where one 
would have Qt instead of da'p (r') as the total TC 
enclosed by dD. However, since Eq. 0 was obtained 
from a continuum generalization, care must be taken in 
cases of continuous p. The discrete nature of TC makes 
p an effective TC density. 

To exemplify the meaning of p in the continuous case 
we produced a radial distribution for simplicity. We con¬ 
sider that p = pqv^^ where po is constant, over a circle of 
radius a and the total TC distributed is Qt- Since now 
p is distributed along a large region, the approximation 
of r ^ 0 is not valid, and the amplitude profile is not 
described by Eq. 0. However, the phase profile from 
Eq. 0 is always valid and equal to 


- 




^ r < a^ 
, r > a. 


( 10 ) 


where we substituted po = Qt {n + 2) /(27ra’^+^) and it 
is assumed that n> —2. 
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Equation (10) is a generalization which smoothly con¬ 
nects usual OV (n = —2) to helico-conical beams, or op¬ 
tical twisters [33l [34] in which n = — 1. Optical twisters 
are interesting because they carry angular momentum 
and also have a higher photon density than the usual 
Laguerre-Gauss or Bessel beams m Therefore they are 
of interest for manipulating particles [34] and may also 
be of interest to nonlinear optics of 0AM carrying beams 
|35| . To our knowledge, other values of n were never pre¬ 
viously reported in the literature. 

We produced beams with the phase profile given by 
Eq. ( p^ , with Qt = 5, fixed a and varying n, and the 
results are shown in Eig. In the phase profiles, Eig. 
[^a-d), it can be seen that larger n increase the phase 
twisting at r < a and reallocates the TC towards the 
border along 0 = tt. The TC displacement can be seen 
also in the zeros of the amplitude profiles, as shown in 
Eig. [^e-h). The LC profiles, Eig. [^i-1), shows that 
larger n values decrease the LC near the center of the 
circle. In Eig. [^m-p) we show the mean LC as a function 
of the radial distance to the center of the circle. The 
values obtained (black dots) agree with the theoretically 
expected from the phase profile (solid blue) via Eq. (10). 
The LC reduction in the center can be understood by 
considering that larger n push p to the boundaries of the 
circle a. Therefore the LC profile is directly associated 
with the p distribution. Also, since the LC is proportional 
to the local classical 0AM of the beam m, this indicates 
that the classical 0AM profile depends similarly on p. 

In summary, we introduced in this work the concept of 
the Topological Potential (TP), Eq. by performing 
conformal transformations over screw dislocations m- 
The identification of the TP paves the way for further 
understanding and tailoring of OV because it creates a 
bridge between OV and 2d electrostatics. Eor applica¬ 
tions where the shape of an OV is relevant, as in optical 
tweezers [8], laser traps [9] or atom guides [10], the TP 
might be used to design OVs for specific applications m- 
Shaped OV may also allow selective excitation of plas- 
monic modes m- However, while the present work can 
be directly used for OV at the focus, further development 
is necessary to understand the effects of propagation and 
address issues as the stability of OV [25]. 

Another important point is that the discussed exam¬ 
ples obtained from the TP might be useful in some appli¬ 
cations. Eor instance, the intensity profile instability of 
TC multipoles may be used to determine the position of 
an extended object image plane, and in aligning spatial 
filters. Another possibility is that, since multipoles form 
a complete (Eourier) basis of orthogonal modes on the 
azimuthal phase, they may be suited for applications in 
quantum communications m- In telecommunications, 
they are an alternative to Laguerre-Gauss beams for data 
multiplexing [12] that can be more stable to turbulence 
[UllT], since the topological information is distributed 
over the beam profile. The 2D TC distributions might be 



n=15 


0 Radius (a.u.) 10 Radius (a.u.) 10 Radius (a.u.) 10 Radius (a.u.) 


Figure 3. Data for 2D radial distribution, Eq. (10), with 
Qt = 5, fixed a and varying n (columns) at z = 0 cm. a is 
graphically represented by cyan dashed lines. The determined 
LC (black dots) at (m-p) corresponds to the azimuthally av¬ 
eraged LC at a given radial distance from the center of the 
circle as a function of the radial distance. The blue solid lines 
in (m-p) correspond to the values expected from the applied 
phase mask. 


used to locally adjust the LC, and consequently the local 
0AM, of a light beam by calculating Eq. @ analytically 
or numerically. This is specially interesting for optical 
tweezers, because then it becomes possible to locally ad¬ 
just the light induced torques. Therefore, one may in 
principle control trapped particles in 2D by simply ad¬ 
justing the TC distribution. 
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